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q ' Abstract. Using high resolution N-body simulations with hydrody- 

^^ . namics and star formation, we investigate interactions and the resulting 

starbursts in galaxies with properties typical of z ~ 3. We apply spectral 
population models to produce mock-HST images, and discuss the ob- 
served magnitude, color, and morphological appearance of our simulated 
galaxies in both the rest-UV and rest-visual bands. 
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^^ ■ 1. Introduction 

^h ; 

OS . An interest in galaxy interactions at high redshift can be motivated from several 

directions. It has been shown that many of the observed properties of galaxies at 

^£3 • z ~ 3-4 (e.g., their number density, luminosity function, colours, sizes, velocity 

Q-f dispersions, etc.) are well-reproduced by a model in which these galaxies are 

predominantly starbursts triggered by galaxy interactions (Somerville, Primack, 
& Faber 1999; SPF). In addition, the collision rate of dark matter halos in high- 
resolution cosmological N-body simulations shows a marked increase at earlier 
epochs, and the clustering properties of these colliding halos at z ~ 3 are similar 
to those of the observed Lyman-break galaxies (Kolatt et al. 1999). 

K> , Previous numerical investigations of starbursts in interacting galaxies (e.g. 

Mihos & Hernquist 1994a; 1996) have assumed initial properties typical of local 
galaxies. Many of these properties (e.g., gas content, surface density, disk-to- 
halo ratio) may be quite different at high redshift. Using the same code as Mihos 
& Hernquist, we are carrying out an ongoing program of simulations to more 
fully explore the parameter space with an emphasis on high-redshift galaxies. A 
more complete version of these preliminary results will be presented in Rosenfeld 
et al. (in prep.). 

2. Simulations 

We use the TREESPH code, including star formation using a Schmidt law (p* = 
Cp^ as ) as described in Mihos & Hernquist (1994b). We set the proportionality 
constant in this equation by requiring the isolated galaxies to lie on the relation 
given by Kennicutt (1998). Masses and mass ratios of typical colliding dark 
matter halos at z = 3 were determined using the simulations discussed in Kolatt 
et al. (1999). For all calculations requiring the assumption of a cosmology, we 
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Figure 1. Star formation rate (left) and observed magnitude at z = 3 
as a function of time since the start of the simulation. In the right panel, 
the bottom curves include an estimate of the effects of dust (see text). 

assume the same cosmological model as those simulations, namely Qo = 0.3, 
i^A = 0.7, and Hq = 70 km/s/Mpc. The mass and exponential scale radius 
of the stars and gas in the galactic disk inhabiting a halo of a given mass at 
the desired redshift are estimated using the semi-analytic models of SPF, and 
are consistent with known properties of Lyman-break galaxies at z ~ 3. The 
disks are assumed to be stable before the start of the interaction. The results 
presented here are for an interaction between equal mass, bulgeless galaxies, 
with halo mass 7.1 x 10 11 Mq, stellar mass 7.2 x 10 9 Mq, stellar scale radius 1.7 
kpc, gas fraction 0.5 and gas scale radius 3.4 kpc. All other properties (relative 
inclination, orbit, etc.), are the same as the fiducial case of Mihos & Hernquist 
(1996; MH96). 



3. Results 

The star formation rate over the course of the interaction is shown in figure 1 
(left panel). The behaviour is qualitatively similar to the results of MH96. We 
convolve this star formation history with stellar population models (GISSEL98, 
Bruzual & Chariot in prep.) to obtain the apparent magnitude of the galaxy 
at z ~ 3 (figure 1, right panel). We have used the solar metallicity models 
with a Salpeter IMF. We show the magnitudes in the WF/PC F606W {V<m) 
and NICMOS F160W (i?i6o) filters (AB system), which probe the rest-frame 
~ 1500A and 4000A part of the spectrum at z = 3. The top set of curves in the 
right panel neglects dust extinction; the bottom set shows the result of including 
a correction of a factor of ~ 5 at 1500 A, as suggested by recent observational 
estimates of typical extinction corrections in bright LBGs (Meurer et al. 1999; 
Steidel et al. 1999), and assuming a Calzetti attenuation curve (Calzetti et al. 
1996). With this level of dust extinction, the bursting galaxy would be visible 
at present spectroscopic limits for about 200 Myr, whereas in the absence of the 
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burst the galaxy would have been well below the detection limit. Note that in 
the absence of dust reddening, the Vejo6 — -^160 color of the galaxy during the 
burst is quite blue (-0.5 - 0), but after the dust correction, typical values (0.5 
to 1.0) are reasonably consistent with observations. 

Figure 2 shows how the interacting galaxy would appear if observed at z ~ 3 
by HST (in the absence of noise or sky background) at various times during the 
merger. Left panels show the WF/PC Vqo6 filter with a pixel size of 0.04 arcsec; 
the right panels show the NICMOS -Hi60 filter with 0.08 arcsec pixels. From top 
to bottom, the galaxies are shown at the beginning of the simulation, during 
the first interaction of the galaxies, after the galaxies have separated again, 
and finally after the final merger. These images can be matched up with the 
star formation and total magnitude curves shown in figure 1 by multiplying the 
simulation time units shown on the figure by a factor of 8.3 to convert to Myr. 

If in fact the observed galaxy population becomes increasingly dominated 
by merging starburst systems at higher redshifts, one would expect to observe 
larger fractions of galaxies with disturbed morphologies and significant substruc- 
ture. We plan to use synthesized images such as figure 2 to develop statistics to 
quantify these effects, and eventually to test whether the morphological proper- 
ties of observed high-redshift galaxies are consistent with the collisional starburst 
scenario. 
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Figure 2. Mock HST images of simulated starburst galaxies at z ~ 3 
at different times during the merger. Left panels shows the Vqoq filter 
(rest ~ 1500A) and right panels show the NICMOS H\qq filter (rest 
~ 4000A). The lightest grey corresponds to a surface brightness of 28 
and the darkest to 18 magnitudes arcsec -2 . 



